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ABSTRACT 

This  paper  describes  the  design  and  use  of  SIGMA,  a  computer  pro- 
gram for  the  calculation  of  radar  cross  section  of  a  dynamic  target 
on  a  pulse  to  pulse  basis,  using  the  MK  25  Fire  Control  Radar.  The 
program  is  written  specifically  for  use  on  the  Scientific  Data  Systems 
9300  computer  in  FORTRAN  IV,  but  is  readily  adapted  to  other  data  pro- 
cessing systems.  The  input  data  base  is  composed  of  magnetic  tape 
recordings  of  instrumented  radar  range  and  signal  strength.  Outputs 
include  target  cross  section  and  range,  and  mean  cross  section  and 
range  in  printed  tabular  listing,  in  graphical  form,  and  on  digital 
magnetic  tape. 
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I.   INTRODUCTION 

Radar  cross  sections  for  various  simple  geometric  objects  have 
been  calculated  and  experimentally  verified.  Cross  sections  for  more 
complex  targets  such  as  fixed  and  rotary  wing  aircraft  have  been 
measured.  This  is  usually  done  by  constructing  a  scale  model  of  the 
target.  The  model  is  then  placed  in  an  anechoic  chamber  and  illuminated 
with  appropriate  RF  energy  at  a  frequency  scaled  to  the  model.  The 
aspect  of  the  model  with  respect  to  the  sensor  is  then  varied  and  the 
backscattered  energy  is  recorded.  From  this  data  cross  section  as  a 
function  of  aspect  and  frequency  is  obtained. 

The  examination  of  cross  sections  of  geometric  objects  is  useful  in 
an  academic  approach  to  diffraction  phenomenon,  but  bears  little  signi- 
ficance to  the  engineering  of  a  radar  system.  Measurement  of  cross 
sections  of  models  is  at  best,  a  first-order  approximation  to  the 
determination  of  back-scatter  cross  section.  The  model  is  static  and 
rests  in  a  controlled  environment.  A  dynamic  target  does  not  show  a 
smooth,  modelled  shape.  Further,  backscatter  is  modified  by  turbine 
engine  modulation  and  changes  in  diffraction  of  the  wave  fronts  due  to 
airfoil  control  surfaces  and  fuselage  vibration.  Since  target  size 
has  a  pronounced  effect  on  the  determination  of  radar  system  performance, 
it  is  then  necessary  that  cross  section  be  measured  on  dynamic  targets. 

One  method  for  dynamic  measurement  of  cross  section  uses  a  ground- 
based  fire  control  radar.  Prior  to  tracking  the  target  to  be  measured, 
the  radar  system  is  calibrated  by  tracking  a  6-inch  metal  sphere  lofted 
by  a  helium-filled  balloon.  Automatic  Gain  Control  (A6C)  voltage  and 


range  are  recorded.  This  calibration  permits  many  of  the  parameters  in 
the  radar  range  equation  to  be  considered  as  a  single  term.  Immediately 
following  the  calibration,  the  target  is  tracked  and  AGC  voltage  and 
range  are  recorded.  Comparing  ranges  to  target  and  sphere  for  equal  AGC 
voltages,  cross  section  for  the  target  can  be  calculated.  While  this 
method  is  a  definite  improvement  over  static  measurement,  disadvantages 
still  exist.  The  AGC  circuit  serves  to  integrate  the  target  data  such 
that  only  mean  cross  sections  can  be  calculated.  Instantaneous  data  is 
lost. 

By  recording  the  pulse  to  pulse  variation  in  signal  strength  and 
computing  a  cross  section  for  each  pulse,  much  statistical  data  can  be 
gained.  For  example,  a  spectral  analysis  of  cross  section  variation 
could  be  performed.  With  a  radar  properly  instrumented  for  recording 
azimuth  and  elevation  data,  correlation  of  the  spectral  lines  with 
pattern  recognition  could  produce  target  signatures.  The  computer  pro- 
gram SIGMA  described  by  this  paper  performs  the  first  of  these  opera- 
tions.  The  received  signal  is  converted  to  radar  cross  section  on  a 
pulse  to  pulse  basis.  The  output  data  is  then  placed  in  a  form  con- 
ducive for  statistical  analysis  by  digital  techniques.  Although  the 
program  was  written  specifically  for  use  on  the  Scientific  Data  Systems 
9300  computer,  it  is  readily  adapted  to  other  user  oriented  systems. 


II.  THEORY 

The  concepts  for  measurement  of  radar  cross  section  on  a  pulse  to 
pulse  basis  are  similar  to  those  outlined  by  Cunningham  [Reference  1]. 
The  fire  control  radar  that  is  used  in  the  measurement  is  first  cali- 
brated by  tracking  a  target  of  known  radar  cross  section,  in  this  case  a 
metallic  sphere.  Since  the  sphere  presents  a  constant  aspect  and  its 
cross  section  may  be  calculated,  a  functional  relationship  between 
radar  range  and  signal  strength  can  be  found.  Armed  with  this  relation- 
ship and  known  constant  sphere  cross  section,  the  unknown  target  is 
tracked,  and  signal  strength  and  range  are  recorded.  Radar  cross  section 
for  the  target  is  ultimately  found  by  forming  a  ratio  of  radar  range 
equations.  Each  range  equation  may  be  reduced  to  a  relatively  simple 
form  if  measurement  of  the  target  immediately  follows  calibration,  and 
atmospheric  effects  are  neglected.  By  proceeding  in  this  manner,  system 
parameters  such  as  transmitted  power,  noise  figure,  integration  improve- 
ment factor,  system  losses  and  antenna  gain,  may  be  assumed  to  remain 
constant.  Let  this  then  be  the  case  as  described  and  further,  let  a 
represent  those  terms  which  remain  constant  between  calibration  and 
measurement.  Then  the  simplified  form  of  the  radar  range  equation  is 
written  as: 

R4  =  ^  (1) 

S 

where:   R  =  radar  range,  meters 

S  =  detected  signal,  watts 

2 
a  =  radar  cross  section,  meter 

o 

a  =  constant,  meter  -watts 


In  the  MK  25  Fire  Control  Radar  System,  the  voltage  output  from  the  box- 
car circuit  is  readily  accessible  and  more  conveniently  measured  than 
actual  signal  power  at  the  receiver  input.  The  functional  relationship 
between  this  voltage  and  signal  power  is  not  linear  because  of  the 
circuitry  involved,  but  may  be  generally  written: 

S  =  f(V)V  (2) 

where  V  =  boxcar  circuit  output  voltage. 
Equation  (1)  may  be  rewritten: 

R4  =  -^  (3) 

f(V)V 


and  further,  let: 


e(v)  =  ^_  (4) 

f(v) 


then  Equation  (1)  becomes 


R4  =  £M£  (5) 


3(V)=^  (6) 


and:  VR? 


s_s 

s 


S(V.)  =  -*-*■  (7) 


where  Equation  (7)  has  been  subscripted  for  the  calibration  data  on  the 
sphere.  Since  Rs  and  Vs  are  measured  and  as  is  calculated,  there  exists 
a  value  of  the  function  which  describes  all  system  parameters.  If  R 
and  V  are  continuously  recorded,  then  the  lumped  parameter  function 
6(V)  becomes  continuous.  For  the  calibration  track,  the  sphere  radar 


cross  section  is  calculated: 

as   =  Tir2  (8) 

where  r  =  radius  of  sphere,  meters. 

It  should  be  noted  that  equation  (8)  applies  only  if  the  circumference 

of  the  sphere  lies  in  the  Optical  region.  The  Optical  region  is  defined 

as  that  region  of  size  where  the  ratio  of  circumference  to  wavelength  is 

equal  to  or  greater  than  10. l 

An  identical  form  of  equation  (7)  is  written  for  the  target  to  be 

measured:  . 

VtRt 

e(v  )  =  J-L.  (9) 

Since  the  same  radar  system  tracks  both  sphere  and  target,  B(VS)  =  $(Vt), 
provided  Vs  =  V  .  Substituting  3(V  )  for  3(Vfc)  in  equation  (9): 

4      4 
R  V   Rt4Vt 

'  =  -5—1  (10) 


a 


S       at 


or:  4 

Rt 

Rs4     S 

To  solve  for  at  it  is  necessary  to  tabulate  Vs  and  R  on  the  cali- 
bration track.  In  a  similar  manner  Vt  versus  Rt  are  recorded.  Since 
B(VS)  must  equal  3(Vt)  for  Vt  =  V  ,  then  there  exists  a  value  of  Rs 


1  See  Skolnik  [Ref.2]  for  a  discussion  of  radar  diffraction  in 
the  Rayleigh,  Resonance  and  Optical  regions. 


for  each  value  of  V.  and  R^.  The  value  of  radar  cross  section  of  the 
sphere  is  known.  Thus,  target  cross  section  may  be  found  directly  for 
each  value  of  V+. 
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III.   INPUT  DATA  BASE 

The  input  data  base  for  both  the  calibration  track  and  target  track 
was  provided  to  the  computer  by  multi-channel  wideband  magnetic  tape 
recording.  Data  extracted  from  the  tape  consisted  of  three  signals: 
boxcar  circuit  output  voltage,  range  marker  ramp  voltage,  and  a  trigger 
pulse.  Since  this  data  must  be  converted  to  digital  form,  a  brief, 
examination  of  the  waveforms  encountered  is  relevant. 

Consider  first  the  dc  level  of  the  boxcar  output.  Output  voltage 
from  the  boxcar  circuit  had  a  wave  shape  for  a  single  pulse  similar  to 
that  shown  in  Figure  (1).  The  ideal (flat)  response  was  corrupted  by 
the  low  pass  characteristics  of  the  tape  recorder  and  its  associated 
matching  networks  at  input  and  output.  It  then  became  necessary  to 
sample  this  signal  at  a  zero-crossing,  delayed  in  time  to  allow 
sufficient  damping.  The  required  delay  was  determined  from  the  recorder 
step  response  and  was  provided  to  the  computer  operator  along  with 
other  parameters  necessary  for  the  analog  to  digital  conversion. 
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Figure  (1) 
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The  radar  range  marker  was  a  ramp  voltage  generated  by  a  servo 
driven  potentiometer  across  a  regulated  dc  power  supply.  The  servo  was 
geared  to  produce  one  complete  revolution  for  each  two  thousand  yards 
range.  The  servo  rotated  clockwise  for  increasing  range  and  counter 
clockwise  for  decreasing  range.  The  range  ramp  voltage  was  biased 
slightly  negative  so  that  the  zero-crossing  could  be  sensed  in  the 
computer  program,  and  used  to  resolve  ambiguities  between  increasing 
and  decreasing  ranges.  The  actual  amount  of  bias  is  not  important  and 
will  be  discussed  in  more  detail  in  Section  IV.  A  typical  wave  form 
for  increasing  range  is  shown  in  Figure  (2).  Figure  (3)  illustrates  a 
decreasing  range.  The  slope  of  the  ramp  is  proportional  to  target 
range  rate. 


Figure  (2) 
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Figure  (3) 
The  trigger  pulse  was  provided  so  that  the  required  delay  could  be 
set  for  sampling  the  boxcar  output.  This  pulse  was  generated  within 
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the  MK  25  Radar  and  was  coincident  with  the  Range  Gate  Enable  signal. 
The  time  between  pulses  was  the  same  as  the  pulse  repetition  period  of 
the  radar,  namely  763  microseconds.  A  typical  trigger  pulse  is  shown 
in  Figure  (4). 
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Figure  (4) 
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IV.  COMPUTER  PROGRAM:  SIGMA 

A.   GENERAL  DEVELOPMENT 

As  described  in  previous  sections,  a  sphere  was  tracked  and  the 
boxcar  circuit  output  voltage,  a  function  of  signal  strength,  and  the 
corresponding  radar  range  were  recorded  for  each  pulse  on  magnetic  tape. 
Similarly,  a  target  of  unknown  radar  cross  section  was  tracked  and  the 
same  measurements,  boxcar  output  and  range,  were  recorded.  It  then 
became  a  matter  of  reading  a  voltage  on  the  target  pulse,  and  since  Vt 
must  equal  V  ,  the  corresponding  range  to  the  sphere  coulti  be  found. 
Equation  (11)  was  then  solved  for  radar  cross  section.  The  most  direct 
scheme  to  perform  these  operations  with  a  digital  computer  would  be  to 
store  in  core  memory  discrete  voltages  and  ranges  for  each  pulse.  In 
effect,  a  table  would  be  constructed  and  a  standard  look-up  routine 
could  be  used  to  extract  the  data.  However,  this  procedure  was  neither 
practical  nor  economical  because  of  the  extreme  size  of  the  input  data 
base.  A  typical  calibration  track,  lasting  fourteen  minutes,  resulted 
in  over  one  million  data  points.  Assuming  that  this  amount  of  data 
could  be  stored,  the  time  required  by  the  look-up  routine  to  extract  one 
point  would  be  prohibitive. 

The  need  for  such  a  cumbersome  table  was  eliminated  by  finding  the 
functional  relationship  between  range  and  boxcar  output  voltage, 
described  in  Equation  (3).  This  was  done  by  calculating  the  mean  signal 
voltage  and  mean  range  for  ewery   1024  pulses,  which  resulted  in  a  new 
data  set  of  approximately  one  thousand  points.  The  functional 
relationship  was  then  found  by  constructing  a  curve  from  the  data  set 
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using  a  least-squared  polynomial  fit,  with  the  optimum  degree  of  fit 
determined  by  application  of  the  F-Ratio  test.  The  polynomial  thus 
found,  together  with  Equation  (11)  were  entered  into  a  single  computa- 
tional program  to  generate  radar  cross  section  on  a  pulse  to  pulse 
basis. 

B.   COMPUTATIONAL  SCHEME 

SIGMA  consists  of  four  separate  programs  written  in  FORTRAN  IV, 
each  containing  subprograms  written  in  either  FORTRAN  IV  or  META-SYMBOL. 
Partitioning  of  the  program  in  this  manner  facilitated  development  in 
that  segments  of  the  main  program  could  be  selected  according  to  the 
amount  of  computer  time  available.  Further,  partitioning  made  possible 
more  rapid  isolation  of  programming  errors  and  shorter  turn  around  time. 
Extensive  use  was  made  of  magnetic  tape  for  input  and  output  operations. 
The  use  of  tape  and  partitioning  was  brought  about  partially  by  the 
limited  core  memory  (32K)  of  the  SDS  9300.  While  appearing  to  be  an 
inconvenience,  the  use  of  tape  allowed  the  results  of  all  data  con- 
versions and  intermediate  calculations  to  be  examined  in  detail.  This 
was  particularly  useful  in  the  analysis  of  noise  spikes  and  their 
subsequent  reduction.  If  desired,  however,  one  can  utilize  the  overlay 
structure  of  the  computer  and  the  four  separate  programs  may  be  com- 
pacted into  a  single  segmented  routine  with  essentially  one  input  tape 
and  one  output  tape. 

1 .  Calibration  Data 

Signals  to  be  processed  are  provided  to  the  computer  on  mag- 
netic tape  in  analog  form  having  the  wave  shapes  shown  in  Section  III. 
In  order  to  reduce  this  data  the  signals  are  first  converted  to  digital 
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form  using  the  COMCOR  Ci  5000  analog  computer,  and  the  Scientific  Data 
System  9300  digital  computer  in  hybrid  configuration.  The  data  is  fed 
into  the  Ci  5000  from  an  Ampex  Model  SP-300  tape  recorder,  sampled  at 
1320  cps.,  digitized  into  binary  words,  and  written  on  magnetic  tape. 
This  entire  operation  is  performed  by  the  program  "A/D". 

Voltages  written  on  tape  by  "A/D"  are  in  the  form  of  24  bit 
word  integers,  scaled  by  amplification  during  the  sampling  process  on 
the  Ci  5000.  It  is  desirable  to  restore  this  data  to  its  original  form 
in  floating  point  notation  for  ease  of  handling  and  output.  Program 
"DATA  SHAPER"  performs  this  function.  In  addition  to  converting  the 
voltages  from  an  integer  representation  to  a  48  bit  real  number,  the 
program  applies  scaling  to  account  for  analog  amplifier  gain.  Further, 
as  the  range  voltages  are  converted,  the  2000  yard  discontinuity  is 
sensed,  and  radar  range  is  generated  directly.  These  ranges  and  their 
corresponding  boxcar  voltages  are  then  stored  and  averaged,  one  record 
at  a  time.  (One  record  contains  1024  data  pairs.)  Both  the  average 
values  and  the  pulse  to  pulse  values  can  be  listed  by  the  line  printer. 
Since  the  record  averages  are  used  to  derive  the  functional  relationship 
between  R  and  V  ,  these  are  written  on  magnetic  tape. 

The  last  phase  of  processing  on  the  calibration  data  is  to  find 
the  polynomial  describing  the  relation  of  R  to  Vs.  Program  "FIT500" 
is  used  for  this  purpose.  FIT500  reads  in  the  record  averages,  selects 
up  to  500  data  pairs,  and  computes  a  least-squares  fit  to  the  degree 
specified  by  the  namelist  input.  Since  the  boxcar  output  is  a  quad- 
ratic function  of  signal  power,  a  second  degree  fit  is  used.  The  pro- 
gram also  executes  an  F-Ratio  test  and  calculates  errors  of  the 
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coefficients  generated.  The  program  outputs  the  polynomial  coefficients, 
the  F-Ratio,  a  tabular  listing  of  the  input  data  with  computed  dependent 
variable  and  error. 

This  concludes  the  processing  of  data  for  the  sphere  track. 
The  coefficients  computed  by  FIT500  are  used  in  the  program  "SIGMAGEN" 
to  calculate  radar  cross  section. 
2.  Target  Data 

Boxcar  circuit  output  voltage  and  radar  range  for  the  target 
being  measured  is  processed  in  much  the  same  manner  as  the  calibration 
data.  The  target  signals  are  first  processed  by  program  "A/D"  for  an 
analog  to  digital  conversion.  This  digital  output  then  serves  as  an 
input  to  the  program  "DATA  SHAPER",  resulting  in  a  magnetic  tape  on 
which  is  written  boxcar  circuit  output  voltage  and  radar  range  in  yards 
for  each  pulse.  This  tape  serves  as  the  input  data  base  for  the 
program  "SIGMAGEN",  which  computes  the  target  radar  cross  section  for 
each  pulse.  SIGMAGEN  provides  a  variety  of  output  options,  including 
both  magnetic  tape  and  line  printer.  In  addition  to  the  pulse  to  pulse 
calculation,  the  program  also  averages  each  record  (1024  pulses),  and 
computes  an  overall  mean  cross  section  for  selected  range  intervals. 
Further,  the  input  tape  need  not  be  started  at  the  load  point,  but  may 
be  advanced  to  any  given  range  by  the  program.  This  permits  small  spans 
of  range  to  be  examined  without  the  need  to  output  an  entire  tracking 
run. 
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C.   IMPLEMENTATION 
1 .  Program  "A/DM 

a.  Analog  Program 

A/D  was  originally  designed  to  perform  an  analog  to  digital 
conversion  on  any  signal  generated  within  the  Ci  5000  with  sampling 
rates  up  to  12,000  sps.  For  the  conversion  necessary  in  computing  radar 
cross  sections,  the  analog  program  shown  in  Figure  (5)  is  used.  Through- 
out this  program  and  all  others,  boxcar  output  voltage  is  labeled  S; 
radar  range,  R;  and  range  gate  trigger,  SYNC.  Both  R  and  S  channels 
from  the  tape  recorder  are  fed  to  separate  50  gain  amplifiers  on  the 
analog  patch  board.  The  gains  are  not  critical  but  should  be  chosen  so 
as  to  provide  a  working  voltage  level  within  the  useful  range  of  the 
Analog  Digital  Converter  unit  (ADC)  in  the  SDS  9300.  The  limits  of  the 
ADC  are  +  99.99  volts.  The  maximum  tape  recorder  output  voltage  is 
approximately  1  volt  RMS.  S  and  R  channel  amplifier  outputs  are 
patched  directly  to  trunks  T500  and  T501  respectively. 

The  SYNC  channel  output  is  patched  through  a  1.0  mfd. 
capacitor  to  a  100  gain  amplifier.  Amplifier  output  is  patched  to  a 
comparator,  the  trigger  level  of  which  is  controlled  by  a  hand  set 
potentiometer.  This  comparator  outputs  a  logic  signal  to  two  delay 
flops.  The  first  delay  flop  is  used  to  set  the  600  microsecond  delay 
of  the  logic  trigger  pulse  with  respect  to  the  SYNC  pulse.  The  second 
delay  flop  is  used  to  adjust  the  width  of  the  logic  trigger  pulse. 
Output  from  the  second  delay  flop  is  fed  to  trunk  T211,  an  interrupt 
line  to  the  digital  computer.  Both  delay  and  width  can  be  set  by 
displaying  the  SYNC  and  logic  trigger  pulses  on  the  CR0.  If  more 
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precision  is  required,  a  period  counter  can  be  used.  To  complete  the 
analog  program,  switch  DS0  is  wired  to  trunk  T100,  a  test  line  to  the 
SDS  9300. 

b.  Digital  Program 

The  digital  portion  of  A/D  consists  of  a  main  program  and 
one  subroutine  written  in  FORTRAN  IV,  and  two  subroutines  written  in 
META-SYMBOL.  Signals  to  be  converted  are  fed  to  the  ADC  on  the  trunks 
specified  above.  If  DS0  is  in  the  ON  position  then  line  T100  has  a 
logic  level  present,  the  interrupts  are  enabled,  and  the  sampling  and 
conversion  process  begins.  Values  of  R  and  S  are  stored  in  array  IBUF, 
dimensioned  (2048).  When  the  array  is  full,  the  values  are  written  on 
magnetic  tape.  This  constitutes  one  record  of  1024  values  of  R  and  S, 
written  in  24  bit  integer  form.  The  actual  value  of  the  integer,  N, 
is  computed  by  the  SDS  9300  according  to  the  following  relation: 

r  .  o23 

N  = —  (12) 

100.0 

Where  E  =  value  of  sampled  voltage  from  ADC. 

The  program  will  continue  to  sample,  convert,  and  output  to  tape  one 
record  at  a  time  until  a  stopping  criteria  is  reached.  To  stop  sampling; 
DS0  can  be  switched  OFF,  or  a  record  limit  can  be  set  and  entered  as 
namelist  input.  A/D  offers  a  number  of  control  options.  The  computer 
will  ask  the  operator  for  an  option  via  the  teletype.  The  option  number 
is  entered  by  the  operator.  The  control  options  are: 

Option  (1)  -  enter  new  namelist  data 

Option  (2)  -  perform  A/D  conversion 

Option  (3)  -  write  "END  OF  FILE"  on  tape 

Option  (4)  -  rewind  tape 
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Option  (5)  -  skip  files 

Option  (6)  -  read  tape,  output  to  line  printer 

Option  (7)  -  read  tape,  output  to  strip  chart  recorder 

The  namelist  inputs,  entered  on  the  teletype,  are: 

NREC  =  total  number  of  records  to  be  written  on  tape 

NCHAN  =  number  of  signal  channels 

NSAMP  =  2048/NCHAN 

ITAPE  =  logical  unit,  magnetic  tape  drive  (integer) 

NDEL  =  number  of  milliseconds  (times  11)  delay  between 
DAC  outputs  for  Option  (7). 

2.  Program  "DATA  SHAPER" 

The  output  from  program  "A/D"  is  in  the  form  of  24  bit  integers 
The  relation  between  the  integer  form  and  the  actual  value  of  the  volt- 
age sampled  is  given  in  Equation  (12).  The  purpose  of  program  "DATA 
SHAPER"  is  to  convert  these  integer  values  to  48  bit  real  numbers  by 
performing  the  inverse  of  the  operations  indicated  in  Equation  (12),  and 
applying  a  scale  factor  to  account  for  gain.  Additionally,  the  program 
computes  the  average  values  for  each  record.  The  input  data  base  con- 
sists of  the  magnetic  tape  output  from  A/D.  This  tape  is  mounted  on 
tape  transport  logical  unit  (1).  Other  input  parameters  are  entered  in 
the  namelist  using  the  teletype.  The  namelist  parameters  are  as 
follows: 

RINT   =  initial  range  to  the  nearest  2000  yard 
multiple,  expressed  as  a  real  number. 

SPANR  =  absolute  value  of  maximum  range  voltage 
minus  minimum  range  voltage. 

SSCALE  =  gain  of  R  channel  amplifier  used  in 
program  "A/D" 


21 


RSCALE  =  gain  of  R  channel  amplifier  used  in 
program  "A/D". 

BIAS   =  absolute  value  of  bias  voltage  on  the 
R  channel  signal . 

FIRST  =  first  record  the  operator  desires  to 
examine. 

LAST   =  last  record  the  operator  desires  to 
examine. 

After  the  namelist  parameters  are  entered,  the  program  will 

read  in  records  from  the  input  tape.  The  records  will  not  be  stored 

until  the  record  specified  by  FIRST  is  reached.  At  this  point,  1024 

values  of  S  and  R  are  placed  in  array  IBUF.  These  values  are  then 

converted  to  48  bit  real  numbers  according  to  the  following  relation: 

E  =   N'100-0  (13) 

223  -XSCALE 

where:      E     =  input  channel  voltage 

N     =24  bit  integer 

XSCALE  =  gain  of  R/S  channel  amplifier 
The  value  of  the  converted  S  channel  voltage  is  stored  in  the  output 
array  BUF.  The  value  of  the  R  channel  is  converted  to  range  in  yards. 
As  seen  in  Figure  (3),  the  servo  driven  potentiometer  causes  a  dis- 
continuity in  the  range  voltage  every  2000  yards.  At  the  discontinuity 
the  voltage  changes  almost  instantaneously  from  a  maximum  to  a  minimum 
for  decreasing  ranges.  Further,  a  small  bias,  0.2  to  0.4  volts,  is 
applied  so  that  the  range  voltage  changes  sign  at  the  discontinuity. 
When  the  program  commences  the  range  is  initialized  to  within  a  2000 
yard  multiple  of  the  range  recorded  at  the  start  of  target  track.  For 
example,  if  the  radar  started  tracking  the  target  at  3600  yards,  the 
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initial  range  (RINT)  would  equal  2000  yards.  If  the  track  started  at 
1650  yards,  RINT  would  equal  zero.  The  computer  adds  or  subtracts  2000 
yards  from  RINT  each  time  the  discontinuity  is  sensed.  The  computer 
determines  the  direction  of  the  range  increment  by  comparing  the  record 
average  just  prior  to  range  voltage  change  from  positive  to  negative,  to 
a  record  average  twenty-five  records  later.  The  actual  value  of  range 
within  the  2000  yard  interval  is  computed  from  the  relation, 


R  =  — 5.   x  2000 


(14) 


where 


Figure  (6) 

VR  =  sampled  range  voltage 

Vg  =  range  bias  voltage 

Vs  =  range  voltage  span. 

Average  boxcar  circuit  voltage  and  average  range  are  computed 
for  each  record  as  it  is  read  from  the  input  data  tape.  These  averages 
are  listed  by  the  line  printer  and  stored  in  core  to  be  written  on  tape 
later,  if  desired.  In  addition,  the  pulse  to  pulse  values  of  voltage 
and  range  can  be  listed  on  the  line  printer,  record  by  record  by 
energizing  Sense  Switch  2.  Pulse  to  pulse  values  can  be  written  on 
magnetic  tape  if  a  tape  is  mounted  on  a  tape  drive  logical  unit  (2). 
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3.  Program  "FIT500" 

The  purpose  of  program  "FIT500"  is  to  derive  the  functional 
relationship  between  radar  range  and  boxcar  output  voltage  for  the 
calibration  track.  The  core  of  the  program  is  subroutine  LSQPL2,  in 
modified  form.  Originally  the  subroutine  was  developed  by  D.E.  Harrison 
in  November,  1969,  to  perform  a  least-squares  fit  using  100  data  points. 
Additionally,  the  subroutine  performed  a  Chi -Square  "goodness  of  fit" 
test  and  an  F-Ratio  test.  In  view  of  the  extreme  size  of  the  input  data 
base  for  program  "SIGMA",  the  subroutine  was  expanded  to  handle  up  to 
500  data  points.  Since  core  storage  in  the  SDS  9300  is  yery   limited, 
LSQPL2  had  to  be  modified  further  by  eliminating  some  of  its  additional 
features.  The  option  to  weight  each  data  point  was  removed,  as  well  as 
the  Chi-Square  test.  The  subroutine  as  it  now  exists  will  fit  500  data 
points  and  output  the  coefficients  of  the  functional  relation.  Other 
outputs  include  a  tabular  listing  of  the  independent  variables,  the 
dependent  variables,  computed  dependent  variables,  and  error. 

Input  data  base  for  FIT500  consists  of  a  magnetic  tape  mounted 
on  tape  drive  logical  unit  (2).  This  tape  should  contain  the  average 
boxcar  voltage  and  average  radar  range  for  each  record  as  generated  by 
the  output  option  in  program  "DATA  SHAPER" .  The  namelist  has  only 
two  entries: 

RECNUM  =  total  number  of  records  on  the  input  tape 
KKM    =  highest  degree  of  fit  desired. 

Normally,  KKM  will  be  set  equal  to  2  whenever  using  program 
"SIGMA",  since  the  relation  between  range  and  voltage  is  quadratic. 
However,  FIT500  can  be  used  alone  to  compute  other  functions  provided 
KKM  does  not  exceed  21 . 
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Calibration  tracks  typically  produced  approximately  1000 
records.  FIT500  reads  up  to  1000  records  and  then  computes  the  function 
using  e^ery   other  data  pair.  In  order  to  read  in  all  records  without 
exceeding  core  storage,  two  subroutine  arrays,  Y  and  DELY,  are  used  in 
addition  to  R  and  S,  the  input  arrays.  Only  the  data  pairs  in  R  and  S 
are  used  for  computation.  The  Y  and  DELY  arrays  are  initialized  to 
zero  before  calling  LSQPL2. 
4.  Program  "SIGMAGEN" 

The  purpose  of  SIGMAGEN  is  to  compute  target  radar  cross  section 
on  a  pulse  to  pulse  basis.  The  input  data  base  consists  of  the  output 
tape  from  program  "DATA  SHAPER",  on  which  are  written  boxcar  output  volt- 
age and  radar  range  from  the  target  tracking  run.  Further,  the  coeffici- 
ents of  the  quadratic  relation  between  boxcar  circuit  voltage  and  range 
for  the  calibration  track  must  be  entered.  As  explained  previously, 
these  are  computed  by  FIT500.  The  input  data  tape  can  be  mounted  on  any 
tape  drive,  the  logical  unit  being  specified  by  INTAPE  in  the  namelist. 

After  namelist  parameters  and  output  options  have  been  speci- 
fied by  the  operator  on  the  teletype,  the  program  causes  the  input  tape 
to  be  read,  one  record  at  a  time.  The  functional  relation  found  by 
FIT500  is  only  valid  over  the  interval  of  ranges  used  in  calibration. 
As  computation  proceeds,  boxcar  voltages  are  read  from  array  BllF,  and 
the  calibration  range,  RCAL,  is  found.  If  RCAL  exceeds  the  calibration 
limits,  it  is  discarded  and  a  new  boxcar  voltage  is  read.  Only  when 
RCAL  is  within  limits  does  a  calculation  for  cross  section  occur.  As 
the  values  for  cross  section,  SIGMA,  are  computed,  they  are  stored  in 
the  output  buffer  array  DBUF.  When  the  array  is  full,  all  1024  values 
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of  SIGMA  and  range  are  averaged.  At  this  point,  output  can  occur 
depending  on  the  option  selected.  There  are  three  different  outputs 
from  this  program.   The  primary  output  of  interest  is  the  pulse  to 
pulse  listing  of  radar  corss  section  and  range.  Secondly,  the  record 
averages  can  be  printed  out.  Finally,  an  overall  average  SIGMA  can 
be  listed.  It  is  conceivable  that  only  a  span  of  specific  ranges  is 
of  interest,  rather  than  the  entire  target  track.  To  provide  for  this, 
the  operator  can  specify  the  starting  range  by  RSTART  in  the  namelist. 
Similarly,  RSTOP  is  entered  as  the  last  range  of  interest.  The  operator 
must  also  enter  the  radar  cross  section  of  the  calibration  sphere  in 
square  meters.  Normally,  a  6-inch  aluminum  sphere  is  used  for  cali- 
bration.  In  this  case  SPHR  should  be  set  equal  to  0.01824. 

It  may  be  desirable  to  examine  the  target  cross  section  even 
though  the  calibrated  range  is  out  of  bounds.  This  can  be  done  by 
executing  Option  (5)  or  (6).  Output  listings  will  then  show  all  cross 
sections  and  ranges,  with  a  minus  sign  before  any  range  computed  out 
of  bounds. 
The  namelist  inputs  are: 

MINLIM  =  range  calibration,  lower  bound 

MAXLIM  =  range  calibration,  upper  bound 

INTAPE  =  input  tape  drive  logical  unit 

OUTAPE  =  output  tape  drive  logical  unit 

2 
SPHR   =  radar  cross  section,  calibration  sphere,  meter 

B1,B2,B3  =  coefficients  of  quadratic  in  ascending  order 

RSTART  =  first  range  in  span  of  interest 

RSTOP  =  last  range  in  span  of  interest 
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The  following  is  a  list  of  options 


Option  (1 
Option  (2 
Option  (3 
Option  (4 
Option  (5 
Option  (6 

Option  (7 
Option  (8 


Output  to  line  printer:  mean  SIGMA  and 
mean  RANGE  for  each  record 

Output  to  line  printer:  SIGMA  and  RANGE 
pulse  by  pulse 

Output  to  tape  and  line  printer:  mean 
SIGMA  and  mean  RANGE  for  each  record 

Output  to  tape  and  line  printer:  SIGMA 
and  RANGE  pulse  by  pulse 

Defeat  calibration  bounds.  Output  to 

line  printer:  SIGMA  and  RANGE  pulse  by  pulse 

Defeat  calibration  bounds.  Output  to  line 
printer:  mean  SIGMA  and  mean  RANGE  for 
each  record 

Read  in  new  namelist 

Stop  the  program 
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V.  CONCLUDING  REMARKS 

Although  SIGMA  provides  the  ground  work  for  computing  radar  cross 
sections  on  a  pulse  to  pulse  basis,  there  are  many  extensions  and  some 
modifications  which  would  expand  the  usefulness  and  facility  of  dynamic 
target  measurement. 

It  is  recognized  that  SIGMA,  as  presently  written  in  four  parts, 
does  not  represent  optimal  programming.  While  partitioning  was  useful 
for  development,  the  program  could  be  consolidated  by  utilizing  seg- 
mented format  in  the  overlay  structure  of  the  SDS  9300.  This  would  also 
reduce  the  need  for  large  amounts  of  digital  tape. 

An  examination  of  the  introduction  of  noise  and  subsequent  error 
generation  would  be  relevent.  The  present  equipment  and  methods  for 
recording  the  analog  tape,  which  is  the  basic  input  to  the  entire 
computer  program,  leave  room  for  improvement.  In  addition  to  tape  hiss, 
wow  and  flutter,  and  switching  transients,  a  certain  amount  of  noise 
is  introduced  by  the  physical  location  of  the  tape  recorder  with 
respect  to  the  radar  equipment  racks.  Some  consideration  should  be 
given  to  the  use  of  a  Dolby  Noise  Reduction  package. 

The  use  of  the  MK  25  Radar  imposes  certain  restrictions.  This 
radar  uses  a  conical  scan  and  therefore  causes  the  returned  signals  to 
be  modulated  at  the  scan  frequency.  While  a  correlation  technique 
could  be  used  in  the  computer  program  to  account  for  scan  modulation,  a 
better  approach  would  be  to  use  a  monopulse  radar.  The  NIKE-AJAX 
radar  system,  for  example,  could  provide  a  better  foundation  on  which 
to  build  an  instrumentation  range. 
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Finally,  regardless  of  the  radar  system  used,  the  computer  program 
could  be  extended  to  handle  target  azimuth,  elevation  and  range  rate. 
Instrumentation  on  airfoil  controls,  as  well  as  target  attitude  should 
also  be  included.  Armed  with  these  inputs,  program  "SIGMA"  could  be 
interfaced  with  an  Adage  Graphics  Terminal  or  similar  display  device 
for  three-dimensional  analysis  of  target  aspect  and  cross  section. 
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